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Abstract 

The critical behavior for intermittency is studied in two coupled one- 
dimensional (ID) maps. We find two fixed maps of an approximate renormal- 
ization operator in the space of coupled maps. Each fixed map has a common 
relavant eigenvaule associated with the scaling of the control parameter of 
the uncoupled one-dimensional map. However, the relevant "coupling eigen- 
value" associated with coupling perturbation varies depending on the fixed 
maps. These renormalization results are also confirmed for a linearly-coupled 
case. 
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A route to chaos via intermittency in the one- dimensional (ID) map is associated with 
a saddle-node bifurcation [|l]. Intermittency just preceding a saddle- node bifurcation to a 
periodic attractor is characterized by the occurrence of intermittent alternations between 
regular behavior and chaotic behavior. Scaling relations for the average duration of regular 
behavior in the presence of noise have been first established by considering a Langevin 
equation describing the map near the intermittency threshold and using Fokker-Plank tech- 
niques. The same scaling results for intermittency have been later found |3J by employing 
the same renormalizat ion-group equation for period doubling with a mere change of 
boundary conditions appropriate to a saddle-node bifurcation. 

Recently, universal scaling results of period doubling for the ID map have been general- 
ized to the coupled ID maps |5|-f|, which are used to simulate spatially extended systems 
with effectively many degrees of freedom [|1(J. It has been found that the critical scaling 
behaviors of period doubling for the coupled ID maps are much richer than those for the un- 
coupled ID map H|[|. These results for the abstract system of the coupled ID maps are also 
confirmed in the real system of the coupled oscillators |TTJ. Similarly, the scaling results of 
the higher period p-tuplings (p —, 3, 4, ...) in the ID map are also generalized to the coupled 
ID maps ||12|| . Here we are interested in another route to chaos via intermittency in coupled 
ID maps. Using a renormalization method, we extend the scaling results of intermittency 
for the ID map to two coupled ID maps. 

Consider a map T consisting of two identical ID maps coupled symmetrically, 

rp . I x n+\ — f\ x n) + g(x n ,y n ), , , 

' \ Vn+l = fiVn) + g{Vm X n ), 

where the subscript n denotes a discrete time, f(x) is a ID map with a quadratic maximum, 
and g(x, y) is a coupling function obeying a condition, 

g(x,x) = for any x. (2) 

The two-coupled ID map (|1]) is called a symmetric map because it has an exchange 
symmetry such that 

a^Taiz) = T(z) for all z, (3) 

where z = (x, y), a is an exchange operator acting on z such that az = (y,x), and a^ 1 is 
its inverse. The set of all fixed points of a forms a synchronization line y — x in the state 
space. It follows from Eq. ([!]) that the exchange operator a commutes with the symmetric 
map T, i.e., oT = To. Thus the synchronization line becomes invariant under T. An orbit 
is called a(n) (in-phase) synchronous orbit if it lies on the invariant synchronization line, 
i.e., it satisfies 

x n = y n for all n. (4) 

Otherwise, it is called an (out-of-phase) asynchronous orbit. 
Let us introduce new coordinates X and Y, 

X = X -V L , Y= X —^L. (5) 
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Then the map (|l|) becomes 



X n+ \ — F(X n , Y n ) 

= l - [f(X n + Y n ) + f(X n -Y n )} 

+ l - [g(X n + Y n , X n - Y n ) + g(X n - Y n , X n + Y n )\, 

(6) 

Yn+i = G(X n , Y n ) 

= l[f(X n + Y n )-f(X n -Y n )\ 

+i [g(X n + Y n , X n - Y n ) - g(X n - Y n , X n + Y n )\. 

This map is invariant under the reflection Y — > — Y, and hence the invariant synchronization 
line becomes Y — 0. Then the synchronous orbit of the old map (P becomes the orbit of this 
new map with Y — 0. Furthermore, the X-coordinate of the synchronous orbit satisfies the 
uncoupled ID map, i.e., X n+ i = f(X n ), because the coupling function g obeys the condition 

(!)• 

Stability of a synchronous orbit of period p is determined from the Jacobian matrix M 
of T p , which is given by the p product of the linearized map DT of the map @ along the 
orbit 

p 



M= l[DT(X n ,0) 



n=l 

" 'f'(Xn) 



n I " f'(X n ) - 2G(X n ) j ' (7) 



n=l 



where f{X) = df(X)/dX and G(X) = dg(X,Y)/dY\ Y=X - The eigenvalues of M, called 
the Floquet (stability) multipliers of the orbit, are 

\i = f[f'(X n ), X 2 = fl[f(X n )-2G(X n )]. (8) 

n n 

Note that Ai is just the Floquet multiplier for the case of the uncoupled ID map and the 
coupling affects only A 2 . 

Consider a synchronous saddle- node bifurcation to a synchronous periodic orbit. The 
synchronous periodic orbit is stable when both Floquet multipliers lie inside the unit circle, 
i.e., |Aj| < 1 for j = 1 and 2. Thus its stable region in the parameter plane is bounded by 
four bifurcation lines, i.e., those curves determined by the equations Xj = ±1 (j = 1,2). 
When a Floquet multiplier Aj increases thorugh 1, the stable synchronous periodic orbit loses 
its stability via saddle-node or pitchfork bifurcation. On the other hand, when a Floquet 
multiplier Aj decreases thorugh —1, it becomes unstable via period-doubling bifurcation. 
(For more details on bifurcations, refer to Ref. [H|.) 



Here we are interested in intermittency just preceding the synchronous saddle- node bi- 
furcation. Employing an approximate renormalization operator P,|T4|- |16[| which includes a 
truncation, we generalize the ID scaling results for intermittency to the case of two coupled 
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ID maps. We thus find two fixed maps of the approximate renormalization operator. They 
have a common relavant eigenvaule associated with the scaling of the control parameter 
of the uncoupled ID map. However, the relevant "coupling eigenvalue" associated with 
coupling perturbation varies depending on the fixed maps. 
Truncating the map (IB) at its quadratic terms, we have 



X n+1 = A + BX n + CXI + FY 2 



[ Yn+i = DY n + EX n Y n 

which is a six-parameter family of coupled maps. Other terms do not appear because 
F(X, Y) and G(X, Y) in Eq. (P) are even and odd in Y, respectively. Here P represents the 
six parameters, i.e., P = (A, B,C, D, E, F). The construction of Eq. (|9|) corresponds to a 
truncation of the infinite dimensional space of coupled maps to a six-dimensional space. The 
parameters A, B, C, D, E, and F can be regarded as the coordinates of the truncated 
space. 

We look for fixed points of the renormalization operator 1Z in the truncated six- 
dimensional space of coupled maps, 

n(T) = AT 2 A-\ (10) 

Here the rescaling operator A is given by 



A= „ . . (ID 




where a is a rescaling factor. 

The operation 1Z in the truncated space can be represented by a transformation of 
parameters, i.e., a map from P = (A, B, C, D, E, F) to P' = (A', B', C, D', E', F'), 

A' = aA(l + B + AC), (12a) 
B' = B{B + 2AC), (12b) 

C = -(B + B 2 + 2AC), (12c) 
a 

D' = D(D + AE), (12d) 

E' = -(BD + D + AE), (12e) 
a 
p 

F'= -(2AC + B + D 2 ). (12f) 

a 

The fixed point P* = (A* , B* , C* , D* , E* , F*) of this map can be determined by solving 
P' = P. We thus find two solutions associated with a saddle-node bifurcation, as will be 
seen below. The map (|^) with a solution P* (Tp*) is the fixed map of the renormalization 
transformation TZ; for brevity Tp* will be denoted as T*. 

For a saddle- node bifurcation at x — 0, the ID map f{x) satisfies 

/(0) = 0, /'(0) = 1. (13) 
Hence the function F(X, Y) in Eq. @ obeys 
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F(0,0)=0, 



dF 
dX 



1. 



(14) 



(0,0) 



We first note that Eqs. ( |12a| )- (|T2c| ) are only for A, B, C, and a. We find one solution 
for A*, B*, C*, and a satisfying the conditions (|14D, 

a = 2, A* = 0, B* = 1, C* : arbitrary number. (15) 

Substituting the values of A*, B* and a into Eqs. ( |12d| )- flT2TD , we have two solutions for 
D*, E* and F\ 



D* = 0, E* = 0, F* = 0, 

D* = 1, £7* : arbitrary number, F* : arbitrary number. 



(16a) 
(16b) 



These two solutions are associated with intermittency in the coupled ID maps, as will be 
seen below. Hereafter we will call each map from the top as the / and E map, respectively, 
as listed in Table |. 

Consider an infinitesimal perturbation e 5P to a fixed point P* of the transformation of 
parameters ( |12a|) -( |l2fl) . Linearizing the transformation at P*, we obtain the equation for 
the evolution of <5P, 



SP' = J5P, 



(17) 



where J is the Jacobian matrix of the transformation at P*. 

Since the 6x6 Jacobian matrix J decomposes into smaller blocks, one can easily obtain 
its eigenvalues. Two of them are 



Ai 



dC 



DC 



A, 



dF' 



OF 



1 + D 



*2 



(18) 



Here Ai is an eigenvalue associated with scale change in X, and hence C* is arbitrary. 
The eigenvalue A2 is also associated with scale change in Y in the case D* = 1; this case 
corresponds to the E map. Thus F* for this case becomes arbitrary. However, in the case 
D* = corresponding to the I map, A2 becomes an irrelevant eigenvalue. Note that the / 
map is invariant under a scale change in Y because F* = 0. 

The remaining four eigenvalues are those of the following 2x2 blocks, 



Mi 



d(A',B') 



Mo 



d(A,B) 
d(D',E' 



d(D,E) 




(19) 
(20) 



The two eigenvalues of Mi (i = 1, 2) are called 5{ and 5' { , and they are listed in Table [Tl| 

The two / and E maps have common eigenvalues of M\. They are 5i = 4 and 5[ = 2, 
which are just the relevant eigenvalues for the case of uncoupled ID maps. Here the 
largest relevant eigenvalue 5i is associated with scaling of the control parameter of the ID 
map near the intermittency threshold. 
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The eigenvalues 62 and 5' 2 of M 2 are associated with coupling perturbations. These 
eigenvalues will be referred to as "coupling eigenvalues" (CE's). The submatrix M 2 for the 
/ map becomes a null matrix, and hence there exist no CE's. On the other hand, the E map 
has a relevant CE 8 2 = 2 and a marginal CE 8' 2 = 1. Here the relevant CE 5 2 is associated 
with scaling of the coupling parameter, while the marginal one 5' 2 is associated with the 
arbitrary constant E*. 

We also obtain the Floquet multipliers A^ and A 2 of the fixed point (0, 0) of the fixed 
map T* of the renormalization transformation TZ. They are given by 

A* = l, \* 2 = D*. (21) 

The I and E maps have a common Floquet multiplier AJ, which is just that for the ID case. 
However, the second Floquet multiplier X 2 affected by coupling depends on the fixed maps; 
A* = (1) for the / (E) map. 

In order to confirm the above renormalization results, we also study the intermittency 
for the linearly-coupled case. The critical set (set of critical points) for the intermittency 
consists of critical line segments. It is found that the I map with no relevant CE's governs 
the critical behavior at interior points of each critical line segment, while the E map with 
one relevant CE S 2 (= 2) governs the critical behavior at both ends. 

We choose f(x) = 1 — ax 2 as the uncoupled ID map in Eq. ([!]) and consider a linear 
coupling case g(x, y) = c(y — x). Here c is a coupling parameter. Three critical line segments 
are found on a synchronous saddle-node bifurcation line a = a c (= 1.75, above which a pair 
of synchronous orbits with period 3 appears. The critical behaviors near the three critical 
line segment are the same. 

As an example, consider a critical line segment including the zero-coupling point c = 
as one end point. Figure [I] shows a phase diagram near this critical line segment denoted 
by a solid line. This diagram is obtained from the calculation of two Lyapunov exponents. 
In case of a synchronous orbit, its Lyapunov exponents are given by 

1 m—1 1 m—1 

a\\{a) = lim - £ ln\f(x n )\, a ± (a,c) = lim - £ In \f'(x n ) - 2c|. (22) 

1 m^oo rn — m^oo rn — 

" L n=0 " L n=0 

Here cry (a±) denotes the mean exponential rate of divergence of nearby orbits along (across) 
the synchronization line y = x. Hereafter, an and a± will be referred to as tangential 
and transversal Lyapunov exponents, respectively. Note also that an is just the Lyapunov 
exponent for the ID case, and the coupling affects only a±. 

The data points on the aj_ = curve are denoted by solid circles in Fig. [I]. A synchronous 
orbit on the synchronization line y = x becomes a synchronous attractor with a± < inside 
the <t_i_ = curve. The type of this synchronization attractor is determined according to 
the sign of a\\. A synchronous period-3 orbit with a\\ < becomes a synchronous periodic 
attractor above the critical line segment, while there exists a synchronous chaotic attractor 
with a\\ > below the critical line segment. These periodic and chaotic regions are denoted 
by P and C in the diagram, respectively. There exists a synchronous period-3 attractor with 
an = on the critical line segment between these two regions. 

The motion on the synchronous chaotic attractor in the region C just below the critical 
line segment is characterized by the occurrence of intermittent alternations between regular 
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behavior and chaotic behavior on the synchronization line. This is just the intermittency 
occurring in the uncoupled ID map, because the motion on the synchronization line is the 
same as that for the uncoupled ID case. Thus, a transition from a regular behavior to an 
intermittent chaotic behavior, which is essentially the same as that for the ID case, occurs 
near the critical line segment joining two end points q = —0.109045 ■ ■ • and c r = on the 
synchronous saddle-node bifurcation line a = a c (— 1.75). 

Consider a "ID-like" intermittent transition to chaos near an interior point with q < 
c < c r of the critical line segment. We fix the value of c at some interior point and vary the 
control parameter e (= a c — a). For e < 0, there exists a synchronous period-3 attractor on 
the synchronization line. However, as e is increased from zero, an intermittent synchronous 
chaotic attractor appears. Like the ID case |§, the scaling relations of the mean duration 
I of regular behavior and the tangential Lyapunov exponent an for an intermittent chaotic 
orbit on the synchronization line are obtained from the leading relavant eigenvalue S\ (= 4) 
of the I map, as will be seen below. 

We first note that the / map is essentially a ID map with zero Jacobian determinant 
(see Table |). Since there exists no relevant CE's associated with coupling perturbation, it 
has only relevant eigenvalues 8\ and 5[ like the ID case. The / map is therefore associated 
with the critical behavior at interior points of the critical line segments. 

A map with non-zero e near a critical interior point is transformed to a new map of the 
same form, but with a new parameter e' under a renormalization transformation 1Z. Here 
the control parameter scales as 

t ' = 5 1 e = 2 2 e. (23) 

Then the mean duration I and the tangential Lyapunov exponent a\\ satisfy the homogeneity 
relations, 

lW) = \m, a ll (e') = 2a ll (e), (24) 
which lead to the scaling relations, 

r(e)~e-", a||(e)~e", (25) 

with exponent 

li = log 2/ log 5i = 0.5. (26) 

The above ID-like intermittent transition to chaos ends at both ends of the critical line 
segment. We fix the value of the control parameter a = a c (= 1.75) and study the critical 
behavior near both ends q and c r by varying the coupling parameter c. Inside the critical 
line segment (cj < c < c r ), a synchronous period-3 attractor with a± < exists on the 
synchronization line, and hence the coupling tends to synchronize the interacting systems. 
However, as the coupling parameter c passes through both ends, the transversal Lyapunov 
exponent aj_ of the synchronous periodic orbit grows continuously from zero, and hence the 
coupling leads to desynchronization of the interacting systems. The synchronous orbit of 
period 3 is therefore no longer an attractor outside the critical line segments, and a new 
asynchronous attractor appears. 
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The critical behaviors near both ends are the same. As an example, consider the case 
of the zero-coupling point c r = 0. Figure || shows the plot of a± versus c for a = a c . Note 
that a± increases linearly with respect to c. Hence a transition from a synchronous to an 
asynchronous state occurs at the zero-coupling end point. 

The scaling relation of er±(c) for a = a c is obtained from the relevant CE 62 (= 2) of 
the E map as follows. Consider a map with non-zero c near the zero-coupling point. It is 
then transformed to a map of the same form, but with a renormalized parameter d under a 
renormalization transformation 1Z. Here the coupling parameter obeys a scaling law, 

d = 6 2 c = 2c. (27) 

Then the transversal Lyapunov exponent a± satisfies the homogeneity relation, 

v ± (d) = 2cx ± (c). (28) 

This leads to the scaling relation, 

v±(c)~c v , (29) 



with exponent 



v 



log2/log5 2 = l. (30) 



Like the case of the I map, the scaling behavior of crii(e) for c = q or c r is obtained 
from the relevant eigenvalue 5i (= 4) of the E map, and hence it also satisfies the scaling 
relation (f25|). The critical behaviors of both exponents cry and a± near an end point are thus 
determined from two relevant eigenvalues 8\ and 62 of the E map. An extended version of 
this work including the results of a renormalization analysis without truncation, the results 
for the many-coupled cases and so on will be given elsewhere [O 
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TABLES 



TABLE I. Fixed point P* of the renormalization transformation 1Z and the rescaling factor a. 



■fivpr] inan 


a A* 


B* 


c* 


£)* p* p* 


± UlCtJJ 


9 n 


i 

_1_ 


Oil Ul Li I OjI V 


n n n 


E map 


2 


1 


arbitrary 


1 arbitrary arbitrary 


TABLE II. 


Some eigenvalues S± 


,s' 1 ,s 2 . 


, and 5' 2 of a fixed 


map T* of the renormalization operator 


are shown. 
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5 2 6' 2 


/ map 
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E map 
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FIGURES 



FIG. 1. Phase diagram of the two-coupled ID map (1) with f(x) = 1 — ax 2 and 
g(x,y) = c(y — x). Here solid circles denote the data points on the a± = curve. The region 
enclosed by the a± = curve is divided into two parts denoted by P and C. A synchronous pe- 
riod-3 (chaotic) attractor with a\\ < (cry > 0) exists in the subregion P (C). The boundary curve 
denoted by a solid line between the P and C regions is just a critical line segment. 

FIG. 2. Plot of the transversal Lyapunov exponent a± of a synchronous period-3 orbit versus 
c for a = a c (= 1.75). 
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